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In the ground state of Ho 2 Ti 2 0 7 spin ice, the disorder of the magnetic moments follows the same rules 
as the proton disorder in water ice. Excitations take the form of magnetic monopoles that interact via a 
magnetic Coulomb interaction. Muon spin rotation has been used to probe the low-temperature magnetic 
behaviour in single crystal Ho 2 _ x Y x Ti 2 0 7 (x — 0, 0.1, 1, 1.6 and 2). At very low temperatures, a linear field 
dependence for the relaxation rate of the muon precession 2(B), that in some previous experiments on 
Dy 2 Ti 2 O y spin ice has been associated with monopole currents, is observed in samples with x — 0, and 0. 1 . A 
signal from the magnetic fields penetrating into the silver sample plate due to the magnetization of the 
crystals is observed for all the samples containing Ho allowing us to study the unusual magnetic dynamics of 
Y doped spin ice. 

I n the spin ice materials R 2 Ti 2 0 7 (R = Ho, Dy) 1-3 a large (~ 10 /.i B ) magnetic moment on the R 3+ ions giving a 
I strong, but at low temperature almost completely screened dipole-dipole interaction, together with a local 
I Ising-like anisotropy leads to an effective nearest-neighbour frustrated ferromagnetic interaction between the 
magnetic moments. The organizing principles of the magnetic ground state in spin ice, or "ice rules", require that 
two R 3+ spins should point in and two out of each elementary tetrahedron in the R 2 Ti 2 0 7 pyrochlore lattice 2,4-7 . 
Excitations above the ground state manifold, which locally violate the ice rules, can be viewed as magnetic 
monopoles of opposite "magnetic charge" connected by Dirac strings 8-10 . Evidence of magnetic monopoles in 
spin ice has recently been observed in several experiments 11-13 . 

Given the existence of magnetic monopoles, it is logical to consider the nature of the magnetic charges and any 
associated currents or "magnetricity". Bramwell et al. used transverse-field muon spin-rotation (TF-/iSR) to 
investigate the magnitude and dynamics of the magnetic charge in Dy 2 Ti 2 0 7 spin ice 14 . In these experiments the 
equivalence of electricity and magnetism proposed in Ref. 8 was assumed and Onsager's theory 15 , which describes 
the nonlinear increase with applied field in the dissociation constant of a weak electrolyte (second Wien effect), 
was applied to the problem of spin ice. It was argued that in spin ice, if the magnetic field B is changed, the 
relaxation of the magnetic moment v„ occurs at the same rate as that of the monopole density and so in the weak 
field limit, v /1 (B)/v /1 (0) = k(B)/k(0) = 1 + b/2, where k is the magnetic conductivity and b = /,i 0 Q 3 B / 8nk 2 B T 2 with 
a magnetic charge Q 14 . At low temperature, the fluctuating local fields lead to a de-phasing of the muon precession 
and an exponential decay in the oscillatory muon polarization as a function of time t 

A(t)=A 0 cos(27tvr)exp(-At), (1) 

where A 0 is the initial muon asymmetry, v = y^BIln is the frequency of the oscillations, and y „ is the gyromagnetic 
ratio. With v / ,(B)/v,,(0) = /l(B)M(0) one can directly infer the magnetic monopole charge. These measurements 
have proven intriguing and controversial. Dunsiger et al.' 6 contend that the TF-^SR data never takes a form where 
k oc v (see however 17 ). It has also been suggested that the magnetic field at any muon implantation site in 
Dy 2 Ti 2 0 7 is likely to take a range of values up to 0.5 T 16 ' 18 ' 19 . If this is the case it is difficult to understand how the 
fields of 1-2 mT used in Ref. 14 could lead to a precession signal. Both Dunsiger et al. 16 and later Blundell 19 have 
suggested that the signals seen in the ^SR data in Ref. 14 originate from outside the sample. In their reply to this 
suggestion, Bramwell et al. 21 acknowledged that their experiments exploited both muons implanted in the sample 
(interior muons) and muons decaying outside the sample (exterior muons), with the aim of separating near and 
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Figure 1 | TF-//SR time spectrum collected at 150 mK in a field of 2 mT 
for a pure Ho 2 Ti 2 0 7 sample. These results are representative of the data 
collected during this study. 

far field contributions to the signal. They went on to note that the 
signal at higher temperatures is dominated by muons implanted in 
the silver backing plate. This possibility was not discussed in their 
original paper 14 . Nevertheless, they continued to insist that the signal 
at low temperature (0.4 > T > 0.07 K) cannot be explained by 
exterior muons and that the Wien effect signal originates from 
muons within the sample or muons sufficiently close to the surface 
of the sample so as to probe the monopolar far field. 

Results 

Fig. 1 shows a TF-/1SR time spectrum collected at 1 50 mK in a field of 
2 mT for a pure Ho 2 Ti 2 0 7 sample. This curve is representative of the 
data collected during this study. A rapid loss in asymmetry from an 
initial value of ~ 0.22 occurs outside the time window of the MuSR 
spectrometer 14-18,19 . The slowly relaxing component of the data were 
fit using Eq. 1. 



Fig. 2 shows the temperature dependence of the muon relaxation 
rate a(T) for Ho 2 - x Y x Ti 2 0 7 extracted from fits to /(SR time data 
collected in 2 mT, (see Methods and Supplementary information). 
For all the samples containing Ho, a nearly T independent 1{T) is 
observed at low-temperature. As the temperature is raised there is a 
rapid increase in 1{T) at some crossover temperature T C r- This T C r 
increases from ~ 0.4 K for the crystals with x = 1.6 and 1.0 (data not 
shown) to 0.5 K for the samples with x = 0.1 and 0.0. Above T CR the 
relaxation rate decreases with increasing temperature and has a sim- 
ilar T dependence for all four samples containing Ho that were 
studied. For two samples (x = 0.1 and 1.6) we also collected field- 
cooled-cooling data. In both cases a divergence between the zero- 
field-cooled warming (ZFCW) and the field-cooled cooling (FCC) 
curves appears at T CR . For pure Y 2 Ti 2 0 7 a temperature independent 
relaxation rate is measured for the whole temperature range (0.05 to 
5 K) studied. 

In order to better understand the origins of these signals we have 
also collected relaxation data as a function of temperature in 2 mT 
for the pure Ho 2 Ti 2 0 7 sample discussed above, covered with a silver 
foil 0.25 mm thick. This thickness of foil is expected to stop all the 
muons before they reach the sample. Muons implanted in silver have 
a negligible relaxation and so any relaxation must result from a 
combination of the externally applied field and/or field lines origin- 
ating from the sample penetrating into the silver. The A(T) curve 
obtained in this way is very similar to the signal from the pure 
Ho 2 Ti 2 0 7 shown in Fig. 2a and demonstrates that at least some of 
the signal comes from fields within the silver, but that these fields are 
the result of the magnetic properties of the sample (see Supple- 
mentary information). 

As a next step we then investigated the magnetic field dependence 
of the muon relaxation rate. Fig. 3 shows X(B) for a sample with x = 0 
at selected temperatures. Studies were also made for samples with x 
= 0.1, 1, 1.6 and 2. Following Bramwell et at, linear fits to the ).{B) 
data were made at each temperature. Using the gradient and inter- 
cept extracted from each fit, the effective magnetic charge Q eff was 
obtained from Q eff = 2.1223 m 1/3 T' 3 , where m = (dl(B)/dB)/A 0 14 . 
For samples with x = 0 and 0. 1 the resulting values of Q e f t range from 
4.5 to 7.5 /ifiA -1 in the temperature regime in which Onsager's 
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Figure 2 | Temperature dependence of the muon relaxation rate A( T) extracted from the fits to the TF-//SR time spectra collected in 2 mT for samples 
of Ho 2 - j:Y x Ti207 with x = 0, 0. 1 , 1 .6 and 2.0. The closed symbols show the zero-field-cooled warming data and the open symbols show the field-cooled 
cooling data. 
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Figure 3 | Magnetic field dependence ofthemuon relaxation rate A(B) for 
pure Ho 2 Ti 2 0 7 at three different temperatures. The values for m = 
(dk(B)l dB)l A a and the effective magnetic charge Q ef j shown in Figure 4 
have been obtained from the straight line fits to the data. 

theory is expected to be valid, but increase rapidly as the tempera- 
tures increase outside this range (see Fig. 4). 

At high temperature, a linear field dependence for a(B) is also 
observed for the two samples with a much higher yttrium doping 
(x = 1 and 1.6) but the calculated Q eff is always greater than 
~ 10 HbAt 1 . For x = 1 and 1.6 in the low- temperature regime T 
< T C r there is no systematic linear field dependence in ).{B) and no 
signal that can be associated with magnetricity. 

We have also looked for a linear magnetic field dependence in a(B) 
for the pure Ho 2 Ti 2 0 7 sample covered in a thick (0.25 mm) silver 
foil. At higher temperatures T > T C r we observed a linear behaviour 
leading to a large Q e ff (i.e. Q e g > 10 [igA^ 1 ), but at low temperatures 
T < T CR we found no signature of magnetricity and could not obtain 
reliable linear fits to the a(B) data or physically acceptable values for 

Qeff 

Discussion 

We can draw a number of important conclusions from our work. Our 
results indicate that at higher temperatures, as suggested prev- 
iously 16,19,21 , the dominant contribution to the a{T) signal arises from 
stray fields from the magnetized spin ice that penetrate into the silver 
sample plate. The observation of a signal in a sample covered with 
thick Ag foil adds weight to this hypothesis. The sample coverage of 
the Ag backing plates used in our experiments was always approxi- 
mately 50%. It will be interesting to explore how this signal changes 
as this coverage is varied. It may also be important to consider the 
ratio between the surface area and the volume of the spin ice in these 
and other experiments. Differences between the bulk and surface 
conductivity of water ice are well documented 20 and it is likely that 
analogous processes operate in spin ice. In reply to the comments on 
their work, however, Bramwell et al. 21 make the point that a signal 
from muons implanted in the sample plate may not negate the 
important findings of their study. Our data are consistent with the 
suggestion made in Ref. 21 that the Wien effect signal may arise from 
inside the sample or from within the Ag sample plate but at distances 
very close to the spin ice sample surface. We will return to this point 
later. First we note that the ).{T) curve for pure Ho 2 Ti 2 07 follows 
closely the form expected for the magnetization of pure spin ice 22 
supporting the view that a(T) reflects the magnetization in all the 
samples studied. This then raises an interesting question concerning 
the low-temperature magnetic dynamics of spin ice. 

Recently there have been a number of experimental reports on the 
magnetic dynamics of spin ice (see for example 23-27 ). In addition to 
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Figure 4 | Q eff versus 1/Tfor samples of Ho 2 _ x Y x Ti 2 07 with x = 0 and 
0.1. The vertical dashed lines indicate the high and low temperature limits 
between which the Onsager theory is expected to be valid 14 and the 
horizontal line marks the value for Q ef j = 4.6 /(bAT 1 8 . The inset shows 
m( T) for the same data; the solid line shows m = Q^ ff / T 2 with 
Q e ff = 5 /i_gA _1 . Also shown in both plots are the data of Bramwell 
et al. from Ref. 14. 

the discussion of magnetic monopoles and the Wien effect 810,14 
authors have also considered the effects of thermal quenching 28 . A 
key component of the current theories of spin ice, is that the magnetic 
response at low temperatures and small applied fields is limited to 
monopole motion. So as the monopole density decreases the char- 
acteristic time scales become longer. This view has recently been 
called into question following new low-temperature AC susceptibil- 
ity measurements that exhibit an activated behaviour with energy 
barriers that are inconsistent with the present understanding of 
monopoles in spin ice 25,29,30 . Our results for the x — 1.6 material, 
showing the survival of ZFCW-FCC splitting in a sample with only 
15% Ho add a further twist to this puzzle. Given the large number of 
non-magnetic "defects" on the corners of many of the tetrahedra in 
this diluted material, it is not easy to attribute the slow relaxation to a 
low monopole density. At such low concentrations of magnetic ions 
even the concepts of a spin ice and monopoles are questionable. 

It is conceivable that single ion physics plays a more important role 
in the behaviour of the diluted materials. Our diffuse neutron-scat- 
tering studies of single-crystal Ho 2 - x Y x Ti 2 0 7 showed that at low 
temperature the scattering patterns are characteristic of a dipolar 
spin ice and appear to be unaffected by Y doping up to at least x = 
1.0 31 . One possible scenario is that effects, such as distortions in the 
local environment due to the variation in the size of the Ho 3+ /Y 3+ 
ions 32 , produce energy barriers at low- T that exceed the cost of an 
isolated monopole. The slow dynamics and the ZFCW-FCC hyster- 
esis at low temperatures would thus cross over from a regime where 
this behaviour is attributed to low monopole density to a regime 
where it is due to exceedingly slow single ion physics. Alterna- 
tively, the long-range nature of the dipolar interactions may give rise 
to collective effects beyond the monopole description which intro- 
duce new energy barriers to spin nipping at very low temperatures 
that occur in both undiluted and diluted systems. The same qualitat- 
ive form for the X(T) data for samples with x = 0.1 and 1.6 indicate 
that additional ingredients may be required to explain the low T 
behaviour in spin ice and that further studies on diluted samples 
are needed to fully understand the role played by factors such as 
impurities, dislocations, and surface effects on the low-temperature 
dynamics of spin ice. 

Returning to the question of magnetricity in spin ice we note that 
in our (iSR data the low-temperature signal that has previously been 
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interpreted as a signature of magnetricity is seen in the x = 0 and 0.1 
samples and is not observed in the more dilute Ho 2 _ x Y x Ti 2 0 7 mate- 
rials. Within the T range indicated by the dashed lines in Fig. 4, where 
the theory presented by Bramwell et al. is expected to be valid, the 
value of Q e ff agrees with expectations. Following Blundell 19 we also 
plot m versus T. We see that the expected m oc T~ 2 only holds for the 
same narrow T range. Our experiments, including two separate runs 
on pure Ho 2 Ti 2 0 7 carried out three months apart, demonstrate the 
reproducibility of the data (see Fig. 2a). A realignment of the 
Ho 2 Ti 2 0 7 disks between runs also shows that the results are not 
particularly sensitive to the exact details of the sample geometry. 
Our results for the samples with a higher Y content and with the 
thick Ag foil demonstrate that the behaviour cannot be attributed to 
instrumental effects. The samples were made at Warwick 33 and are 
Ho rather than Dy based pyrochlores, eliminating the possibility of 
material specific results. 

In summary transverse-field ;<SR experiments on Ho 2 - x Y x Ti 2 0 7 , 
including measurements on non-magnetic Y 2 Ti 2 0 7 and a sample of 
Ho 2 Ti 2 0 7 covered in thick silver foil, suggest that the majority signal 
in the A(T) response comes from stray fields due the sample mag- 
netization penetrating into the silver sample plate 1619 . The results for 
Ho 2 Ti 2 0 7 are comparable with those observed for Dy 2 Ti 2 0 7 . The 
low- temperature (T < T CR ) linear field dependence in /1(B) is only 
observed in samples with x = 0 and 0.1. In this low-temperature 
regime the value of Q eff agrees quantitatively with the theory pre- 
sented in Ref. 14. The low-temperature hysteresis in 1{T) for the 
magnetically dilute material (x = 1.6) appears inconsistent with 
the current understanding of monopoles in spin ice. 

Methods 

Single crystals of Ho 2 _ x Y x Ti 2 0 7 (x = 0, 0.1, 1, 1.6 and 2) were grown in an image 
furnace using the floating zone technique 33 . The single crystal disks were glued on to a 
silver plate and covered with a thin (0.01 mm) sheet of silver foil to improve thermal 
conductivity. The plate was then attached to the cold stage of an Oxford Instruments 
3 He/ 4 He dilution refrigerator. Transverse-field muon spin-rotation experiments were 
performed using the MuSR spectrometer at the ISIS pulsed muon facility, Rutherford 
Appleton Laboratory, UK. The magnetic field was applied along the [001] direction, 
perpendicular to the initial direction of the muon spin polarization which was along 
the [110] axis. Measurements were carried out as a function of applied field at fixed 
temperature and as a function of temperature in a fixed magnetic field. See 
Supplementary information for full details of the measurement protocols. 
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